Acoustic radiation controls dynamic friction: 
Evidence from a spring-block experiment 
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Brittle failures of materials and earthquakes generate acoustic/seismic waves which lead to radi- 
ation damping feedbacks that should be introduced in the dynamical equations of crack motion. 
We present direct experimental evidence of the importance of this feedback on the acoustic noise 
spectrum of well-controlled spring-block sliding experiments performed on a variety of smooth sur- 
faces. The full noise spectrum is quantitatively explained by a simple noisy harmonic oscillator 
equation with a radiation damping force proportional to the derivative of the acceleration, added to 
a standard viscous term. 



The science of solid friction has a long history, dat- 
ing back in the western world to the geometrical work of 
Leonardo di Vinci, continuing with the empirical Amon- 
tons' law two centuries later and Coulomb's investiga- 
tions of the influence of sliding velocity on friction in the 
18'th century. Some 40 years ago, R. Feynmann in his 
famous lectures stated that "Friction is a very compli- 
cated matter... and in view of all the work that has been 
done it is surprising that more understanding of this phe- 
nomenon has not come about." 

Only three decades ago was it recognised that friction 
plays a role in the mechanics of earthquakes and it 
was proposed that stick-slip as observed in friction || 
was relevant to earthquake dynamics. Numerous lab- 
oratory experiments have been carried out to identify 
the relevant parameters controlling the solid friction and 
hence the stick-slip behaviour ||. Low velocity (below 
« 1 cm/s) experiments have established that solid fric- 
tion is a function of both the velocity of sliding and of 
one or several state parameters, characterising the true 
surface of contact Q|. These so-called Ruina-Dieterich 
laws now constitute the basic ingredients in most mod- 
els and numerical elasto-dynamic calculations directed at 
understanding the process of rupture and earthquake nu- 
cleation. 

A well-known and serious limitation of these calcula- 
tions based on laboratory friction experiments is that the 
friction laws used have been determined under steady- 
state sliding conditions using velocities of no more than 
« 1 cm/s (and often much less), i.e., considerably below 
the sliding velocity of m/s occuring during an earthquake. 
Thus, one may ask whether it is correct to extrapolate 
these laws and their velocity weakening dependence to 
higher velocities relevant to earthquakes? Such consider- 
ations become all the more relevant when one examines 
the underlying physical mechanisms of the friction laws. 



At low velocity, effects such as hysteretic elastic and plas- 
tic deformations at the scale of roughness asperities seem 
to play a dominant role j5j. At larger velocities, new 
mechanisms come into play. Collisions between asper- 
ities and transfer of momentum between the directions 
parallel and perpendicular to the motion are potentially 
important mechanisms ||. Recently, Tsutsumi and Shi- 
mamoto Q have performed friction measurements on ro- 
tating cylindrical samples at velocities up to 1.8 m/s and 
for slips of several tens of meters. Their results indicate 
a change of regime from velocity weakening to velocity 
strengthening at large velocities. This is confirmed by 3-d 
numerical simulations performed in a regime of velocities 
of meters to tens of meters per second Q , reflecting the 
increasing strength of vibrational damping. 

Radiation damping is well-documented in electromag- 
netism [^,0 and nuclear physics. Surprisingly, the corre- 
sponding mechanism of dynamic friction due to radiation 
of phonons or seismic waves has received little attention 
(see however JTT| ) , notwithstanding its large potential im- 
pact on the dynamics. For instance, it is well-known that 
Burridge-Knopoff spring-block models does not recover 
the correct elasto-dynamic continuous limit but rather 
lead to a Klein-Gordon equation with a mass term im- 
plying finite range interactions [jOJ. This problem can 
be addressed by adding a viscous damping accounting 
for radiation losses with an amplitude finely tuned to 
the critical damping value JT~3[ | . The influence of ultra- 
sounds on crack dynamics in brittle materials has been 
demonstrated by using both the natural sound emitted 
by the propagating crack and an artificially generated ul- 
trasound burst |L4[ . Although the acoustic energy is only 
5% of the energy needed to propagate the crack, the pres- 
ence of sound waves in the specimen strongly modifies the 
fracture dynamics because the sound interacts with the 
crack tip. 
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Here, we re-analyse the high-frequency part of the 
power spectrum from an experimental investigation of 
stick-slip in dry metallic friction We focus on the 

dynamic friction in the high velocity regime, i.e., slip 
velocities ranging up to v max sw 0.35 m/s, which pro- 
vides a direct demonstration of the role and nature of 
radiation damping feedback on the dynamics. Other as- 
pects of these experiments have been reported elsewhere 
Jl5| , [l6[ , but a satisfying explanation for the high velocity 
behaviour could not be proposed at that time. 




FIG. 1. Top view of the experimental setup. A small 
steel block (1) is placed on a rotating steel table (2) and via 
a thin plastic rod attached to a steel spring. The steel spring 
is placed inside a metal box (3) kept fixed in the laboratory 
frame (4). The metal box (3) forms a common shield with 
the electronic measuring devices and the battery (not shown) 
powering the Wheatstone bridge containing the strain-gauges 
thus reducing electro-magnetic noise. Details of the setup 
contained in the metal box (3) is shown in figure EJ. 




FIG. 2. Side view of the metal box containing the spring 
and the strain gauges glued on each side of the spring. 

Details of the experimental setup can be found in 
|l5l , fL6| and here we will just repeat the essential features, 
see figures [l] and ||[ The experiment consisted of a small 
steel block (1) placed on a rotating steel table (2) and 
via a thin plastic rod attached to a steel spring inside 
a metal box (3), see figure [l] and caption. Other metal 
surfaces were used and gave similar results. The other 



end of the steel spring was kept fixed in the laboratory 
frame of reference (4) . As the surface of the table moved, 
the spring became elongated until the spring force ex- 
ceeded the friction force on the block and the spring con- 
tracted until the friction force again exceeded the spring 
force and so forth. The mass of the block was 28.5 g 
and the spring constant 22.4 N/m giving a natural fre- 
quency of the spring and block ljq = yjk/m k, 28 s~b 
The elongation of the spring was measured using strain 
gauges mounted on each side of the spring, see figure 
^, and placed in a Wheatstone bridge. The use of bat- 
teries both for the bridge and the pre-amplifier as well 
a common shield reduced electromagnetic noise signifi- 
cantly and the over-all signal-to-noise ration was better 
than 1000:1, corresponding to an elongation of the spring 
of 5/im to 5 mm or more. 

When the block moves relatively to the surface of the 
table, there are many collisions between the asperities of 
the two surfaces. Now, the asperities typically have a 
spacing of 10 — 100 microns, which means that the block 
(on average) will experience thousands of collisions be- 
tween asperities per second. The motion of the block can 
thus be modeled as a "noisy" damped harmonic oscilla- 
tor with the following equation of motion for the position 
x(t) of the block 

x + c rad x +c vis x + ujI'X = ri(t), (1) 

where rj(t) is a "white" noise term, i.e., (f](t)) = and 
(Tj(t)rj(t')) = bS(t — t') accounting for the stochastic mo- 
tion of the asperities. We introduce two damping terms. 
The non-standard c ra d x term is the Abraham-Lorenz 
expression for the first-approximation of the direct reac- 
tion force due to radiational damping. A viscous friction 
c V i S x must also be present in order to produce a friction 
force at constant velocities. It includes the effect of all 
other friction mechanisms, including a renormalisation of 
the Abraham-Lorenz term at the scale of individual as- 
perities that undergo acceleration/deceleration even un- 
der constant block velocity. As we shall see below, the 
radiative term dominates completely above the charac- 
teristic frequency loq. 

It is useful to recall the derivation of the radiation 
damping reactive force F rac i — c rac i x . We start from 
the general expression for the power P(t) radiated by 
this element. By Galilean invariance, P(t) must be zero 
if the velocity x is constant and becomes non-zero when 
the acceleration x is non-zero. Assuming analyticity and 
symmetry under x — > — x and performing a Taylor ex- 
pansion in powers of x, we get the leading term as 

P(t) = m e Tx 2 . (2) 

The quadratic dependence of the radiated power by a 
small acceleration element is so general that it applies 
to any physical problem involving an accelerating body 
coupled to a wave. m e is the mass of the element and r is 
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a characteristic time proportional to R/c where R is the 
typical linear size of the element and c is a wave velocity. 
Expression (||) recovers the Larmor power formula for 
the electromagnetic radiation of an electric charge. It 
also describes the acoustic radiation from an accelerating 
volume element, the fluid gravity waves radiated from an 
accelerating surface distortion or even the gravitational 
waves from an accelerating black hole. 

In order to obtain the expression of the reaction force 
F ra ,d due to radiation, we follow Jackson B and view, 
by the requirement of energy conservation, the radiated 
power P(t) as minus the work per unit time of F rac [ : 
dt F ra d x = — dt P(t). Integrating the r.h.s. by 
part and neglecting the boundary term (which are zero 
for periodic motion), we get the Abraham-Lorenz expres- 
sion for the radiative force 



Frad — TTl e T X 



(3) 



As a first-order approximation, it must be replaced by 
an integro-differential equation when radiation damping 
becomes the dominant term in the dynamics. 

The consequence of this result is dramatic. For a given 
oscillatory amplitude, the radiative damping is propor- 
tional to u! 3 compared to the usual to for viscous damping. 
This corresponds to a weaker damping at low frequency 
and a more efficient effect at large frequencies (but still 
sufficiently small so that the wavelength remains larger 
than the source size). 

In (|l]) , a renormalised coefficient c ra d is used instead of 
m e r jm to account for the collective behaviour summed 
over all asperities. We can estimate it theoretically by 
transforming expression (|2|) in the Fourier domain as 
Pu = c ra d,muj 2 V 2 , where m is the block mass and V 
the particle velocity of the sound wave generated by the 
moving mass. We then equate this expression to the ra- 
diated power of N coherent asperities of radius R and 
contact pressure p given by 1 .2N 2 oj 2 (nR 2 p) 2 / pc 2 [p7| , 
where p « 7.8g/cm 3 , c « 5800m/s for steel and the co- 
efficient 1.2 is for Poisson ratio of 1/4 and varies slowly. 
Notice that NirR 2 p ps mg, assuming that the contacts 
are represented by discs of radius R, where g — 9.8 m/s 2 
is the earth acceleration. This approximation assumes 
that there is no significant additional inertial pressure 
due to vertical motion of the block other than its static 
weight on the contact asperities. This leads to the fol- 
lowing simple expression : 



1.2 / mg 

pm \ Vc 



(4) 



All parameters in (Q) are known except the particle wave 
velocity V, which can be determined from our fit to the 
spectrum. 

The power spectrum corresponding to (|l|) is 

2b 

S ( U1 ) = —^ ^72 —, 72 > ( 5 ) 



leading to an ~ 1 /w 6 -decay of the power spectrum for fre- 
quencies larger than the natural frequency luq w 28 s _1 . 
Figure || shows the power spectrum measured experimen- 
tally in the regime where the block was constantly sliding. 
Note the corner- frequency (or shoulder) at ujq ps 28 s _1 . 
There appears to be a 1 ju) 2 background below this corner 
frequency. This is most likely an artifact of making a fi- 
nite time measurement thus creating an illusion of a slow 
constant drift in the signal. Indeed, if x(t) — > x(t) + at 
then x(u>) — ► x(u) + a/iuj (for periodic boundary condi- 
tions) resulting in the addition of a a 2 /o; 2 -term to equa- 
tion (|). 

In figure |§| is also shown a fit to the data with equation 
(H) with this correction term added and with only viscous 
damping, i.e., (c ra d — 0). It emphasises the importance 
of the radiation term to account correctly for high fre- 
quency tail of the spectrum. In figure || we also show a 
fit of the high frequency tail of the spectrum with the pre- 
dicted l/w 6 -decay from the radiation term. The agree- 
ment between the data and the prediction is excellent. 
The fit of the entire spectrum with both damping terms 
non-zero in equation (^|) is found to be highly unstable 
due to the complete dominance of the radiative damping 
compared to the viscous damping indicating that most 
of the useful high-frequency spectrum shown in figure || 
is controlled by the novel radiation damping term. 



Eq. (5) without radiation - 
'Measured power spectrum' - 
l/omega"6 - 




(u 2 ~ W 2 ) +UJ 2 (c r adU 2 



Angular frequency 

FIG. 3. The power spectrum of the motion of the block 
in the experimental setup shown in figures ^ and |j| The fit 
is equation (^) with a a/w 2 -term added to account for the 
apparent drift and with c ra d = and uo = 28s ~ . The high 
frequency tail of the spectrum has been fitted with a pure 
w -6 -decay. 



From the fit of the high frequency tail of the spectrum 
we get c r ad ~ 0.02 s. Using (|^) we obtain an estimation 
V ~ 25 (Ltm/s for the wave particle velocity giving an 
acoustic pressure pcV w 1100 Pa corresponding to very 
energetic local sources. 
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This finding may illuminate the so-called "cut-off" fre- 
quency problem on the source acceleration spectra 
of moderate to strong Californian earthquakes. The un- 
expected result is that the seismic spectrum falls off very 
fast beyond a frequency attributed to a characteristic size 
of cohesive fault size or to scale-length of heterogeneities 
of the fault plane. Our results may quantify this phe- 
nomenon by confirming the role of the radiation by accel- 
erating asperities and specifies the quantitative shape of 
the spectral fall-off. However, earthquakes faults have a 
thick layer of fault gouge between the sliding rock blocks 
which might affect the magnitude of the radiation term. 
In order to clarify this issue, further experimental work 
with such an intermediate layer should be performed. 

In conclusion, we have presented a new analysis of 
high sensitivity measurements of noise spectra in spring- 
block experiments |Te|| with sliding velocities in the range 
« 0.35 m/s. We have shown that the spectrum ex- 
hibit an approximate l/o; 6 -decay for frequencies larger 
than the natural corner value ljq = yjk/m. This in- 
dependence can be rationalised simply by the generic 
Abraham-Lorenz radiation damping law, with a reac- 
tive force proportional to x , generated by the radiation 
of sound waves due to collisions of asperities. We ex- 
pect that this finding extends to the slip-stick regime 
for which measurements of the two-point correlations of 
successive slip characteristics, slip distance and time of 
slip were found to be very weak |l6| ] indicating highly 
nonlinear dynamics. The implication of this finding for 
rupture and earthquake modelling is of great potential 
impact. Our results suggest that the Ruina-Dieterich 
friction laws cannot be extended to the high velocity 
regime relevant for earthquakes and many cases of rup- 
ture. Future works include the generalisation of the 
acoustic Abraham-Lorenz radiation damping beyond the 
first-order approximation and the derivation of a gener- 
alised Fluctuation-Dissipation theorem relating the noise 
amplitude b to the damping coefficients c ra d and c V i S . 
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